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Abstract 



The problem of the experimental determination of the amount of entangle- 
ment of a bipartite pure state is addressed. We show that measuring a single 
observable does not suffice to determine the entanglement of a given unknown 
pure state of two particles. Possible minimal local measuring strategies are 
discussed and a comparison is made on the basis of their best achievable 
precision. 



PACS-numbers: 03.67.-a, 03.65.Bz 



Typeset using REVTeX 



I. INTRODUCTION 

Quantum mechanical states of multiparticle systems can be entangled, a fact well known 
since the early days of the formulation of the quantum theory [p. However, the status of 
this property has changed substantially in recent years. Entanglement can also be viewed 
as a resource and as such features in different processes of potential practical importance, 
for example quantum teleportation f2j, quantum cryptography || or even high precision 
measurements |4]] . From a theoretical point of view, entanglement of bipartite pure states, 
and its properties under local quantum operations, are reasonably well understood. There is 
a unique measure of entanglement for these systems, provided by the von Neumann entropy 
H, and optimal ways for entanglement manipulation are known |J. However, there is a 
remaining practical question which has not been addressed so far: How can one optimally 
measure the amount of entanglement of an unknown bipartite pure state? 
At a first sight, this question may seem obvious. Reconstructing the reduced density operator 
of any of the two subsystems will do the job. However, the essential point is that we require 
the determination to be optimal and the reconstruction of the reduced density matrix may 
provide redundant information, given that we are asking for just a feature of the composite 
state, its entanglement. This is a single number and the first question to be answered is 
whether there exists a single operator whose experimental measure may provide us with just 
the amount of entanglement of the state. Note that further details of the state itself are not 
of interest in the problem we are posing here J7J. We will prove in the following that such 
an operator does not exist, a conclusion that confirms what it could initially be thought as 
an educated guess. Knowing the impossibility of a test using a repeated measurements of 
a single observable, we will discuss possible strategies aimed to be minimal, in the sense of 
involving the smallest number of observables. In order to avoid any ambiguity when counting 
the number of observables involved in a given measurement protocol, we will define such 
number as the different number of meters each observer has to read out. Among different 
minimal strategies, that is, strategies involving the same number of meters, we will call 
optimal the one providing the best accuracy when supplied with the same resources. We 
will show that, in fact, measuring the reduced density matrix turns out to be a minimal way 
to proceed. Moreover, the protocol can be made optimal when involving projections along 
mutually orthogonal directions. We have organized the paper as follows. In section II we 
state formally the problem. Section III shows the impossibility of finding a single observable 
whose measure may allow the experimental determination of the amount of entanglement. 
Minimal local strategies are discussed in Section IV. When supplied with the same number 
of identically prepared bipartite pure states, we discuss the performance of two classes of 
minimal measurements from the point of view of the achievable precision in determining the 
amount of entanglement. Section V is devoted to conclusions. 

II. AN EXPERIMENTAL SCENARIO 

Let us imagine the following situation. We are provided with a state preparator which 
creates pairs of two-level particles (qubits) in an unknown entangled state. These entangled 
pairs are distributed to two remote locations where two observers, Alice and Bob, may 



perform local measurements as well as interchange classical communication. The internal 
dynamics of the device is not specified and the only thing Alice and Bob know is that, with 
high accuracy, the state they share is pure. Therefore, the two-qubit state can be written as 

P=\i/>ab)(iI>ab\, (1) 

where 

I ^)ab = a o I 00) + a x | 01) + a 2 | 10) + a 3 | 11) . (2) 

In this expresion (| 0) , | 1)) refer to the eigenvectors of the operators a z , the complex co- 
efficients di, (i — 0, ... ,3), being completely unknown. In addition, we assume that the 
machine may supply a large number of identical pairs. The aim is to use the resulting pairs 
for a quantum information task and, therefore, the only property we are interested in is its 
amount of entanglement. Moreover, we require the measurement aimed to determine the 
amount of entanglement to be optimal in the following sense. First, the protocol should 
involve the smallest possible number of observables. Such tests will be called minimal. And 
secondly, among minimal tests, we will define as optimal the class of protocols that yield the 
best resolution when supplied with the same resources, i.e., the same number of identically 
prepared two-level systems. 

The problem is still rather general and, for simplicity, three further assumptions will be 
made: 



1. The experimental situation is such that it only allows to act on one pair at a time. 
In other words, we restrict ourselves to incoherent measurements. Alice and Bob are 
not allowed to store a given number of particles and perform a joint measurement on 
them @. 

2. No ancillary systems are available and the only allowed incoherent measurements are 
projective ones. 

3. The adopted protocol is rigid, in the sense that we will not accumulate information 
from a given set of initial measurements and re-adjust our strategy afterwards. 

In these conditions, we will show that no single operator measurement allows to determine 
the amount of entanglement of an unknown bipartite pure state. 

III. IMPOSSIBILITY OF A SINGLE-OBSERVABLE MEASURING STRATEGY 

The amount of entanglement of a bipartite pure state is given by its von Neumann 
entropy, 

E(iPab) = -tr(p A log 2 p A ) = -tr(p B log 2 p B ) (3) 

where Pa(b) = t r B(A)P is the reduced density matrix of each subsystem and p is given by 
Eq. (1). In terms of the concurrence C 0, defined as 
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the amount of entanglement can be expressed in a compact form as follows 

i + Vi-c 2 \ 



(4) 



£(V> 



ABj 



1 + vT 1 


-c 2 


2 


1- Vi- 


-C 2 



log 2 



log 2 



(5) 



1- vi -c r 



It should be noted that if all coefficients a^ were real, the concurrence could be obtain via 
the repeated measurement a single observable, a y <8> cr y . We will now prove that in general, 
i.e. where no a priori information is provided about the state of the bipartite system, it 
is not possible to evaluate C 2 (i/jab) by means of measuring a set of orthogonal projectors 
JPi — I Oi) ( Oi\, J2i= IPi = 11, where the | Oj)'s form an orthonormal basis of certain operator 
O. This measurement would allow us to compute the four probabilities pi = \(Oi \ ip) | 2 and 
therefore it provides three independent real numbers. It is obvious that this will not be 
enough to fully reconstruct the pure state | ip) but one may still ask whether the resulting 
information may be enough to compute a property of the state, its amount of entanglement. 
In order to check this, let us first re-write the concurrence C in a more convenient form. For 
that, we will express the state | iPab) in terms of the eigenbasis of the operator O as 
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where the coefficients m^'s are purely real and ipi G [0, 2tt). Then, C 2 can be written as 
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Let us define a new matrix K with elements given by 



K i:j = (Oi\a y ®a y Of] 



In terms of this quantity, the squared concurrence can be written as 



(7) 
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Looking at this expression one can already formulate the guess that it will not be possible 
to obtain C 2 from just measuring the probabilities Pi, given that no information about 
the relative phases ipi will be unveiled by the measurement. In what follows we will prove 
explicitly that C 2 , and therefore the amount of entanglement of the bipartite pure state, 
cannot be an univaluated function of the probabilities pi (equivalently, of the coefficients 



rrii 



A. An useful lemma 

Let us define two new auxiliary matrices S and a with matrix elements given by 

S ij = {0* j \O i ) (10) 

and 

Vij = (Oi\(Xy ®dy \Oj) . (11) 

It is easy to check that the following properties hold. 

1. The matrix K of Eq. 8 satisfies K = K T , as it follows immediately from the hermiticity 
of the operator a y ® a y . 

2. If the | Oj)'s form an orthonormal basis, the corresponding conjugate vectors | O*) also 
form an orthonormal basis. Then, the matrix S defined above is just the change of 
basis matrix between the two representations, i. e., 

\O l ) = T,^\0*), (12) 

3 

and therefore S*S = 1 and \Det(S)\ = 1. In particular, this implies that Det(S') ^ 0. 

3. Det(a) = 1. 

We now have all the ingredients for proving the following lemma. 

Lemma: det(K) ^ . 

Proof: The matrix elements of K can be written in terms of those of o and S as follows: 

O* 
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Therefore, K = aS^ and the lemma follows from properties 2 and 3. 

B. Impossibility of a single-observable test 

We will now prove that assuming that the measurement of a single observable allows to 
determine the concurrence of the state, and therefore its amount of entanglement, yields 
a contradiction with the previous lemma. Given that the state | ip) is unknown, the test 
we are seeking must be universal, that is, the hypothetical observable O has to provide the 
amount of entanglement of whatever input state. The idea underlying our proof is to show 



that there will always be a particular case yielding to a contradiction. Therefore, if no a 
priori information is provided, the minimal test will necessarily require measuring more than 
one observable. 

Consider the particular case where m = Vfi\ = l/v2 and m 2 = m 3 = 0. In this case, Eq. || 
takes the form: 
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are linearly independent yields the set of equalities 
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Eq. [15] implies that either Kqo = and/or K\\ = 0. Taking Eq. ^ into account, this 
corresponds to the cases where K i ^ or K 01 = 0. In other words, we obtain that two 
among the three complex numbers (K 00 , K 01 , K u ) must be zero. 

If we repeat this argument for all the cases where any two of the coefficients rrii are equal 
to 1/a/2 and the remaining two equal to zero, we end up with the requirement that in all 
the following sets of three complex numbers 



(K 00 ,K 0U K U ) 
(K m K 12 ,K 22 ) 
(K 2 2, K23, K 33 ) 



(-^oo? K 02 , K 22 ) 
(K U ,K 13 ,K 33 ) 



(K 00 ,K 03 ,K : 
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there must be at least two of them equal to zero in any set. This fact imposes a certain 
symmetry for the allowed i^-matrices. Explicitly, K can only be one of the following 
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and analogous forms obtained when interchanging the roles of the indexes, or of the form 



K 3 



( K Q1 K 02 K 03 \ 

K Q1 K 12 K 13 

K 02 K 12 K 23 

\ K 03 K 13 K 23 / 



(18) 



It should be noted that many other cases could be obtained if any of the matrix elements 
written as non-zero were zero, however these additional cases are not of interest here, as will 
become clear afterwards. 

Our proof ends with showing that in any of the allowed forms for K, some of the possibly 
nonzero coefficients in K turn out to be zero. Therefore, all the allowed forms for K, i.e. all 
forms compatible with the requirement of being the concurrence a univaluated function of 
the real numbers pi, will have determinant equal zero, which contradicts the lemma stated 
before. 

This can be easily shown for matrices of the form K\ or K 2 just following an argument 
parallel to one used above and choosing three of the coefficients m; equal to l/v3 and the 
remaining one equal to zero. Let us analyze here the case of i^-matrices of the form K 3 . 
If we set the rrij's coefficients to the values m = and rrii = l/\/3 for (i = 1, ...,3), the 
squared concurrence given by Eq. reads 
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where we have introduced the relative phases a = <p 3 — <p 2 , (3 = <p 3 — (j)% and 7 = <p 2 — <f>\. 
Using again the argument invoked in proving the allowed forms for the matrix K, we obtain 
that the following equalities must hold 



K 12 K* 3 = 
K 12 K* 3 = 

K 13 K* 3 = 0. 



Therefore, within the three numbers (Ki 2 ,Ki 3 ,K 23 ), two must be zero. If we repeat the 
argument making each time one of the m^'s zero and the other three equal to a/3, we end 
up with allowed forms for K 3 with either four or two matrix elements different from zero. 
For instance, an allowed form of K is given by 
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( K m \ 

K 01 

K 23 

V K 23 / 



(19) 



But choosing now all the coefficients m, equal to 1/2 suffices to get the constraint ^01-^23 = 
0, which yields a zero determinant for K. The same reasoning applies to the other five 



possible cases. Therefore, assuming that C 2 is only a function of the probabilities pi yields to 
the condition det(K) = 0, but we have proved in the previous section that the orthonormality 
of the vectors Oj demands the determinant of K to be non zero. As a result, it is not 
possible to find a single operator O whose measurement allows to determine the amount of 
entanglement of the pure state | iPab)- 

IV. MINIMAL TESTS 

The previous analysis shows that a measuring strategy employing a single-observable 
does not allow Alice and Bob to know the amount of entanglement of the state they are 
sharing. Knowing this fact, the natural question to ask is to determine the minimal mea- 
suring strategy than may allow them to evaluate E(iPab)- It is clear that if they measure 
two different observables O ab of the form analyzed before, they will fully reconstruct the 
original pure state and can, therefore, compute its amount of entanglement. It is a remark- 
able fact that acting on the whole Hilbert space of the two particles, we cannot isolate the 
information related to the amount of entanglement alone by means of measuring a single 
observable. If no information is known about the state, determining its amount of entangle- 
ment leads to a full reconstruction of the state. However, such a non-local implementation 
may not be the easiest to implement experimentally and local strategies are preferred. We 
will discuss in this section possible ways to proceed if one is constrained to act locally and 
we will determine the expected precision of the protocols. 

A. Local actions without exchange of classical communication 

When Alice and Bob are constrained to act locally and no classical communication can 
be exchanged among them, the minimal measuring strategy corresponds to the local recon- 
struction of the reduced density operators. For instance, Alice may reconstruct the operator 
Pa = t r BPAB and send at the end of the protocol a final message to Bob whose length will 
depend on the required resolution. 

To achieve this, Alice needs to perform three projective measurements along linearly inde- 
pendent directions. We will show in the following that measuring along three orthogonal 
directions is in fact optimal, in the sense that choosing this configuration yields the smallest 
associated uncertainty in the experimental determination of the determinant of the reduced 
density operator (See Eq. (4)). 
Let us write the reduced density matrix in the general form 

^ i+ ^=\u + + %, s i-ty (2o » 

in terms of the corresponding Bloch's vector. With the above parametrization we have a one 
to one correspondence between directions in three dimensional space and directions within 
the Bloch sphere. Note that the determinant of pa only depends on the modulus of Bloch's 
vector. In other words, it is rotationally invariant. Suppose now we are planning to measure 
the amount of entanglement projecting the reduced density matrix of a given state along 



three linearly independent directions. The uncertainty associated with this measurement 
will depend on: 

1. the modulus of the corresponding Bloch's vector, as the amount of entanglement does. 

2. the relative position of Bloch's vector with respect to the three projective directions. 

Because of conditions 1 and 2, assuming the initial distribution of states to be isotropic, the 
average uncertainty after measuring sufficiently many states will only depend on the relative 
position of the three directions we project along. In particular, this implies we can choose a 
given direction to be the z-axis. We will call h and rh the other two directions so that the 
angles they form with the z-axis are 6 n and 9 m . Then we can write the average uncertainty 
as 

®av J \"ni "mi H'nra) • 

Here, <p nm is the relative azimut angle <p m — <f) n . Moreover, because of condition 2, the 
following equalities must hold 

J \"ni "mi H'nra) J \7T "ni "mi " ^Pnm) i 

f(0ni ®mi 4>nm) = f(Pni ^ ~ &mi <P 
f{0rn Qm, 4>nm) = /(tT ~ 
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This is equivalent to say that we could redefine the three positive axis (simultaneously or 
not) without changing the average uncertainty. Finally, from the previous set of equations 
one obtains that the function / has to have an extremum at 



Or, 



vr/2, n = 71-/2, nm = 7r/2. 



We have numerically calculated the average uncertainty 5 av over 10.000 states of the compos- 
ite system uniformly distributed over the complex four dimensional joint Hilbert space and 
calculated this uncertainty for a series of different measurements with 6 m = ir/2, 6 n = n/2 
but <f) nm ranging from to n. Results are shown in Figure 1, where we show the mean 
uncertainty in determining the amount of entanglement S av , defined as 
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(21) 



as a function of the relative difference nm for a fixed value of m = n . In the definition 
above, Pq, P m and P n are the probabilities to obtain a spin up when measuring along 
directions z, rh and h respectively and where each SP denotes the squared variance given 
by 5P = - N '- . The bracket means the average over the isotropic distribution of states. It 
is clear that 5 av reaches in fact a minimum when <p nm = ir/2. Similar figures can be plotted, 
all of them supporting that the minimum uncertainty is indeed achieved when the three 
directions of projection are chosen to be mutually orthogonal. 
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FIG. 1. Average uncertainty as a function of the relative azimut angle (p mn for a fixed value of 
&m = dm- The average uncertainty in determining the amount of entanglement is minimal when 
the three linearly independent directions of projection are chosen to be orthogonal. 



B. Local actions with exchange of classical communication 



Let us now assume that Alice and Bob agree to cooperate. Then, the amount of entan- 
glement can be evaluated from the measurement of two Pauli operators in each side. If they 
agree to measure different operators in each round, they again fully reconstruct the state. 
However, if one of them always measures the same Pauli operator, for instance they choose 
to compute the observables a z <8> I and 1 <g> a z and, in a subsequent round, the observables 
era;® 1 and l<g>a z , for which they should read out three meters, they can obtain the amount of 
entanglement but they will neither get full information about the state itself, nor about the 
whole reduced density matrix. Indeed, if we denote by Pi (i=0, 1,2,3) the four probabilities 

associated to the outcomes (++,H — , — h and ) when measuring a z <g> 1 and \®a z and by 

P ++ ,P + ^,P^ + and P the corresponding probabilities when measuring a x <g> 1 and 1® <r z , 

it can be easily shown that the probabilities Pij (i,j = +, — ), can be written in terms of the 
probabilities Pi as follows: 



10 



P++ = ^\P + Pi + 2^P P X cos0 ol ) (22) 

P+- = \ (Po + Pi ~ 2/^ cos 0oi 

1 



P-+ = g ( 1 -Po-Pl + 2\/P 3 P 2 COS 23 
1 



P— = -[l-P -Pi- 2^/ P 3 P 2 COS (f) 23 

where we have rewritten the amplitudes of the initial state as a* = rriie 1 ^', (i=0,l,2,3), and 
called (fiij = 4>i—(j)j. From the previous set of equations, we see that the functions cos(0 o — 0i) 
and cos(02 — 03) can be expressed in terms of measurable quantities in the form 

cos(0 o - 0i) = — ++ /TT - fT - (23) 
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Noting that 



C 2 = 4 P X P 2 + P P 3 - 2JP P 1 P 2 P 3 cos(0 o - 0i + 03 - <h) • (24) 



we see that this measuring strategy suffices to determine the squared concurrence and there- 
fore the amount of entanglement of the pure state but it does not allow the full reconstruction 
of the initial state. As we will show in the next section, this protocol is not optimal, in the 
sense of providing the best accuracy when measuring locally a minimal set of observables. 



C. Which strategy yields the best resolution? 

The measuring strategies described above are both minimal, in the sense of involving the 
smallest number of meters to be read out. However, it is not obvious whether the precision 
achieved following these two strategies will be the same. In fact, we will show in the follow- 
ing that it is not. When provided with the same resources, that is, using the same number 
of identically prepared entangled pairs, we can get the amount of entanglement with higher 
precision by means of the local reconstruction of the reduced density operator. If we denote 
by iV the number of entangled pairs, iV being large in the statistical sense, a numerical 
simulation with 10 6 states from an isotropic initial distribution yields the following results. 



1. The measurement procedure by means of the local reconstruction of the reduced den- 
sity operator has an associated uncertainty which scales with iV as 

0.3 

Oloc ~ 



N 
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2. The associated uncertainty with a local measurement of the form described in Sec. 



IV.B is substantially much larger ||10|| . More precisely 

2.3 



Oloc+cc ~ ^- 

Note that, once N is given, the resulting number of measurements in each measuring protocol 



is different. While in the first case each single probability will scale as P ~ l/JN/3, the 
larger number of measurements yields each probability in the second procedure to scale as 
Pnl/jN/2. 

From these results one may be lead to the conclusion that the best resolution will always be 
achieved by means of reconstructing the reduced density operator. However, this may not 
be true. Imagine that, in the context of the second protocol, Bob measures a different Pauli 
operator. If the direction of projection is orthogonal to the z-axis, this procedure will also 
allow to reconstruct the initial state. Will now the associated uncertainty be reduced with 
respect to the case analyzed above? In the light of the results obtained when measuring 
the reduced density operator, this should be the case. It should be noted, however, that 
the number of observables required in this measurement protocol is no longer minimal, as it 
would require the observer Bob to read out an additional meter. 

V. CONCLUSIONS 

We have analyzed the problem of determining experimentally the amount of entangle- 
ment of bipartite pure states when one has a large supply of identically prepared systems 
on which one is restricted to act by means of projective measurements. We showed that, 
provided that the entangled state is totally unknown, no measuring strategy involving a 
single operator exists. Therefore, acting on the Hilbert space of the composite system does 
not allow to single out the amount of entanglement without allowing to determine the state 
completely. When local actions are considered, the minimal protocol for determining the 
amount of entanglement involves measuring three different observables. We have analyzed 
here two classes of minimal tests. In the first one, no exchange of classical information is 
required and entails the local reconstruction of the reduced density operator. The procedure 
is optimal, in the sense of having the smallest associated uncertainty, when measuring along 
three mutually orthogonal directions. The second class of protocols requires the use of clas- 
sical information. Here we have analyzed a possible strategy and showed that it suffices to 
determine the amount of entanglement of the pure state but not its full reconstruction. The 
associated resolution turns out to be worse than the one corresponding to the measurement 
of the reduced density operator. The analyzed protocol is not necessarily the most precise 
among the whole class of measuring strategies by means of local actions with the exchange 
of classical information. However, an increase in the resolution would be done at the prize 
of increasing the number of meters to be read out and the protocol would no longer be 
minimal. Establishing how the best accuracy achievable within a protocol of this type, and 
which allows the full reconstruction of the state, compares with the one associated to the 
reconstruction of the the reduced density operator is an interesting open problem. 
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